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Abstract. Stabilization is still a somewhat controversial issue concerning its very 
existence and also the precise conditions for its occurrence. The key quantity to set- 
tle these questions is the ionization probability, for which hitherto no computational 
method exists which is entirely agreed upon. It is therefore very useful to provide var- 
ious consistency criteria which have to be satisfied by this quantity, whose discussion 
is the main objective of this contribution. We show how the scaling behaviour of the 
space leads to a symmetry in the ionization probability, which can be exploited in the 
mentioned sense. Furthermore, we discuss how upper and lower bounds may be used for 
the same purpose. Rather than concentrating on particular analytical expressions we 
obtained elsewhere for these bounds, we focus in our discussion on the general princi- 
ples of this method. We illustrate the precise working of this procedure, its advantages, 
shortcomings and range of applicability. We show that besides constraining possible 
values for the ionization probability these bounds, like the scaling behaviour, also lead 
to definite statements concerning the physical outcome. The pulse shape properties 
which have to be satitisfied for the existence of asymptotical stabilization is the van- 
ishing of the total classical momentum transfer and the total classical displacement 
and not smoothly switched on and off pulses. Alternatively we support our results by 
general considerations in the Gordon- Volkov perturbation theory and explicit studies 
of various pulse shapes and potentials including in particular the Coulomb- and the 
delta potential. 



INTRODUCTION 

There is considerable interest in the high intensity regime (intensities larger 
than 3.5xWcm -2 for typical frequencies), because since the early nineties it may be 
realized experimentally. The perturbative description, which was a very successful 
approach in the low intensity regime, breaks down for such high intensities. Thus, 
this regime constitutes a new challenge to theorists. Comparing the status of 
the understanding and clarity of the description of the two regimes one certainly 
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observes a clear mismatch and should probably conclude that the challenge has not 
been entirely met so far. One also observes a clear imbalance between numerical 
calculations and analytical descriptions. 

In particular, the issue of stabilization has led to several controversies and there 
are still several recent computations which are in clear contradiction to each other. 
Since it is not very constructive simply to count the numbers of numerical results 
which agree and those which do not 1 , our investigations aim at analytical descrip- 
tions which unravel the physical assumptions and might serve to pinpoint possible 
errors. 

In view of the panel discussion at this meeting the main purpose of this con- 
tribution is to summarize our findings [1-6] and in particular explain the working 
and limitations of our method in the hope to dispel a few misunderstandings and 
misconceptions which have occurred. 



FRAMEWORK AND PHYSICAL PROPERTIES 

We start by stating our physical assumptions. We consider an atom with poten- 
tial V (x) in the presence of a sufficiently intense laser field, such that it may be 
described in the non-relativistic regime by the time- dependent Schrodinger equa- 
tion in the dipole approximation 

i dj^ = ^_* +v{ fl + s .fi ( ^ ^ t)=H (f> t) ^ t y (1) 

We will use atomic units throughout this article. We take the pulse to be of the 
general form 

E(t) = E f(t) (2) 

where f(t) is assumed to be a function whose integral over t is well behaved with 
f(t) = unless < t < r. This means r constitutes the pulse duration, f(t) 
the pulse shape function and E the amplitude of the pulse, which we take to be 
positive without loss of generality. 

Important quantities for our discussion are the total classical momentum transfer 
b(r), the classical displacement c(r) and the classical energy transfer a(r) defined 
through the relations 

b(t) = [ dsE{s), c(t)= [ dsb(s), a{t) = - [ dsb 2 (s) . (3) 
Jo Jo 2 Jo 

The quantity of interest, which one aims to compute, is the ionization probability 
V ((f) defined as 



' Panel discussion at this meeting. 



V M = \\(1-P)U (r, 0)<p\\* = l- \\PU (r, 0)<p\\\ (4) 

Here P denotes the orthogonal projection in the space L 2 (R 3 ) of square integrable 
wave functions onto the subspace spanned by the bound states if of H (x, t — 0) , 
||-|| is the usual Hilbert space norm and the time evolution operator is defined by 

U (t, t'^j = T[Exp(-i J H{x, s)ds)} , (5) 

with T denoting the time ordering. The question one is interested in is: How does 
V(ip) behave as a function of E 7 In particular is it possible that V(<p) decreases 
when the field amplitude E increases, in other words does stabilization exist? 
Quantitatively this means we should find a behaviour of the form 

^Mi^)<0 for V M *1 (6) 

with < Eq < oo on a finite interval for Eq. We refer to a behaviour in (6) for 
the equal sign as weak stabilization and for strict inequality we call this strong 
stabilization. 

We stress once more that this description is entirely non-relativistic. The rel- 
ativists regime surely poses a new challenge and a full quantum field theoretical 
treatment is desirable, but it should be possible to settle the question just raised 
within the framework outlined above, since stabilization is not claimed to be a 
relativistic effect. In particular it is not clear which consequences on the physics 
in this regime one expects from a description in the form of the Klein-Gordon 
equation 2 . Furthermore, appealing to a more formal description 3 in terms of scat- 
tering matrices 4 instead of the time evolution operator U (t,t) will not shed any 
new light on the question raised, unless one deals with non-trivial asymptotics. 

The time-ordering in (5) poses the main obstacle for the explicit computations of 
V ((p). To get a handle on the issue, one can first resort to general arguments which 
provide analytical expressions constraining the outcome. The least such arguments 
are good for is to serve as consistency checks for results obtained by other means. 
This is especially useful when one has a controversy as in the case at hand. In 
addition we will demonstrate that they also allow some definite statements and 
explain several types of physical behaviour without knowing the exact expression 
of the quantities which describe them. 



CONSTRAINTS FROM SCALING PROPERTIES 

More details concerning the arguments of this section may be found in [5]. Denot- 
ing by A > the dilatation factor and by 77 the scaling dimension of the eigenfunc- 

2 ) See contribution to the panel discussion at this meeting by F.H.M. Faisal. 

3 ) See contributions to the panel discussion at this meeting by F.H.M. Faisal and H. Rciss. 

4 ) For pulses of the form (2) the scattering matrix S = limt ± ^± QO cxp(it + H + ) ■ U(t+,t-) ■ 
cxp(-it-H-) and U(t, 0) coincide in the weak sense, (see e.g. [1] for a more detailed discussion) 



tion (p(x) := ip(x, t — 0) of the Hamiltonian H (x, t — 0), we consider the following 
scale transformations 5 

x — > a? = Xx and (/9(f) — > ^'(a 7 ') — A _, V(:r) . (7) 

As the only two physical assumptions we now demand that the Hilbert space norm, 
i.e. = ||</?'(x , )||, remains invariant and that the scaling of the wavefunction 

is preserved for all times. From the first assumption we deduce immediately that 
the scaling dimension has to be 77 = d/2 with d being the dimension of the space. 

The scaling behaviour (7) may usually be realized by scaling the coupling con- 
stant. Considering for instance the wavefunction <p(x) = \fa exp(— a\x\) of the 
only bound state when the potential in (1) is taken to be the one-dimensional 
delta-potential V(x) = aS(x), equation (7) imposes that the coupling constant has 
to scale as a — > a' = X^a. Choosing instead the Coulomb potential in the form 
V(x) = a/r requires the same scaling behaviour of the coupling constant for (7) to 
be valid. This is exhibited directly by the explicit expressions of the corresponding 
wavefunctions (p n im(x) ~ a 3 ^ 2 (ar) 1 ex.p(—ar/n)L^t^(2ar/n). 

From the second assumption we conclude 

ijj(x , t) -> V t') = U'(t', 0)<// (£' ) = A- d/ V(f , t) = X- d/2 U(t, O)ip(x) . (8) 

Consequently this means that the time evolution operator should be an invariant 
quantity under these transformations 

U(t u t ) = T (e-'ti H ^ ds ) - U'(t' 1: Q = T (V* J S> H ' { "' s)d ^ = U(t u t ) . (9) 

Equation (9) then suggests that the scaling of the time has to be compensated by 
the scaling of the Hamiltonian in order to achieve invariance. Scaling therefore the 
time as 

t -> t' = xn , (10) 

equation (9) only holds if the Stark Hamiltonian of equation (1) scales as 

H (x, t) -> H' (x', t') = X m H (x, t) with r} H = -r} t . (11) 

The properties (10) and (11) could also be obtained by demanding the invariance of 
the Schrodinger equation (1). The overall scaling behaviour of H (x, t) is governed 
by the scaling of the Laplacian, such that we obtain the further constraint 

VH = -2 . (12) 

5 ) More formally we could also carry out all our computations by using unitary dilatation opera- 
tors U(X), such that the transformation of the eigenfunction is described by U(\)ip(x) = \ n ip'{\x) 
and operators O acting on ip(x) transform as U(X)OU(X)^ 1 = O' . 



As a consequence we can read off the scaling properties of the potential as 



V (x) -> V (f ) = \~ 2 V (x) 



(13) 



Considering for instance the one-dimensional delta-potential and the Coulomb po- 
tential in the forms specified above, equation (13) imposes that the coupling con- 
stant has to scale as a — > a' = \~ 1 a in both cases. This behaviour of the coupling 
constant is in agreement with our earlier observations for the corresponding wave- 
functions. 

We will now discuss the constraint resulting from equation (11) on the scaling 
behaviour of the pulse. We directly observe that 



This equation is not quite as restrictive as the one for the potential, since in the 
latter case we could determine the behaviour of the coupling whereas now a certain 
ambiguity remains in the sense that we can only deduce 



E ^E' = \^oE , /(t) -/'(*')= with VE +Vf = -3- (15) 



Thus, under the assumptions we have made, it is not possible to disentangle the 
contribution coming from the scaling of the amplitude or the pulse shape function. 
However, there might be pulse shape functions for which rjf has to be 0, since 
no suitable parameter, analogously to the coupling constant for the potential, is 
available in its explicit form to achieve the expected scaling. 

Finally, we come to the scaling behaviour of the ionization probability. Noting 
that the projection operator has to be a scale invariant quantity, i.e. P — > P' = P, 
we obtain together with (7) and (9) that the ionization probability remains an 
invariant quantity under the scaling transformation 



We have therefore established that transforming the length scale corresponds to 
a symmetry in the ionization probability V (tp). This symmetry can be exploited 
as a consistency check in various approximation methods in numerical or analytical 
form as outlined in [5]. In this sense the arguments of this section are similar in 
spirit to those of the next section. Nonetheless, scaling properties may also be used 
to explain directly certain types of physical behaviour, as for instance the behaviour 
of V (<f) as a function of the coupling constant (see [5]). 



E (t) -> E' (t) = \~ 3 E (t) 



(14) 



V M = \\(1-P)U (r, 0) tp\\ 2 -> V {ipf) =V(<p). 



(16) 



CONSTRAINTS FROM BOUNDS 



In this section we wish to comment on the method of computing bounds which 
is alternative to computing V(<p) exactly. This means we estimate expressions of 
the form 



\\(1-P)U(t,0M 2 <V u &) and ||Pf/(r, 0)y?|| 2 < l-Vi(<p) (17) 
such that 

VM < r(<p) < ?M ■ (is) 

How does this work? We can not go into all the technical details, but we would 
like to illustrate the general principle of the computational steps involved. First 
one should note that from a mathematical point of view there are seldom general 
principles for deriving such inequalities, except for a few elementary theorems (see 
e.g. [7]). Therefore the steps in the derivations very often do not always appear 
entirely compelling. In mathematics, absolute inequalities, i.e. those which hold for 
all real numbers, are important in analysis especially in connection with techniques 
to prove convergence or error estimates, and in physics they have turned out to be 
extremely powerful for instance in proving the stability of matter [8] or to establish 
properties of the entropy [9]. 

The basic ingredients which are always exploited are the Minkowski and Holder 
inequalities 

W + VW < W + OT, < UW-Wl (19) 

used in the form 

U -^ + X-X\\< U -x\\ + \\X - ^'||, (20) 

Uxx-^'wkuxw-wx-^'w , (21) 

where ip and ip' are meant to be formal objects. The aim and sometimes the art of 
all considerations is now to choose X such that the loss in accuracy is minimized. 
One should resort here to as much physical inspiration as possible, for instance if 
there is a conjecture or a result from other sources which suggests a dynamics one 
can compare with. There exist also more sophisticated possibilities to estimate the 
norm, as for instance to relate the Hilbert space norm to different types of norms, 
e.g. the operator norm 6 or the Hilbert-Schmidt norm 7 

ll^ll < \\A\\ op U\\ < U\\h.s.U\\ ■ (22) 

Where do we start? In fact, the starting point is identical to the one of perturba- 
tion theory, that is the Du Hamel formula involving the time evolution operator 
associated to two different Hamiltonians Hi(t) and H 2 (i) 

U 1 (t,t') = U 2 (t,t')-i f dsU l (t,s)(H 1 (s)-H 2 (s))U 2 (s,t) . (23) 

Jt' 

6 ) The operator norm is denned as |j^4|| op =sup ¥>: || ¥> || = i 

7 ) Denoting by ot\ > 012 > ■ ■ ■ the positive eigenvalues of the operator T = (A* A) 1 / 2 the Hilbert- 
Schmidt norm of the operator A is defined as ||-A||#.s. = E^Li a «) 1,/2 - 



For instance, identifying the Stark Hamiltonian in (1) with Hi(s), one chooses 
H 2 (s) = —A/2 + x ■ E (t) or H 2 (s) = -A/2 + V (x) in the high- or low intensity 
regime, respectively. Instead of iterating (23) and ending up with a power series in 

V in the former or a power series in E in the second case one inserts (23) into (17) 
and commences with the estimation of the norm in the way just outlined. Most 
conveniently these considerations are carried out in a different gauge, for the high 
intensity regime in the Kramers-Henneberger gauge. 

Where do we stop? The whole procedure may be terminated when one arrives at 
expressions which may be computed explicitly. 

When can we apply bounds? In general in all circumstances. In particular problems 
occurring in the context of perturbative considerations, like the convergence, are 
avoided completely. Especially when the strength of the potential and the field are 
comparable, e.g. in the turn-on and off region, this method is not limited in its 
applicability, as is for instance the case for the Gordon- Volkov series. 
What can we deduce? Ideally Pi(<f) and P u (<f) are very close to each other, in 
which case we are in the position of someone solving the problem numerically 
with Pi(<p) and P u (y?) related to the numerical errors. If the lower bound tends 
to 1 for an increasing finite realistic value of Eq there will be little room left for 

V ((f) to decrease and one may deduce that stabilization is absent (see figure 9 in 
[2]). Furthermore, we can always make statements about the extreme limits. For 
instance for the extreme frequency limit we obtain 



This relates our discussion to the seminal paper on the stabilization issue by Gavrila 
and Kaminski [10]. For the extreme field amplitude limit we found 



where ipGv{ T ) — Ugv{ t -, 0)y? is the Gordon- Volkov wave function. For the definition 
of Ugv see (39). We would like to stress that this limit is not merely of mathematical 
interest 8 . The result (25) is a clear indication of weak stabilization, though it is 
still desirable to find the precise onset of this behaviour. It is also clear that as a 
consequence of (25) a value of V (tp) which is equal or larger than the r.h.s. of (25) 
for any finite and experimentally realisable value of Eq immediately implies the 
existence of strong stabilization. 

What are the shortcomings? For realistic values of the parameters involved the 
expressions sometimes yield 




(24) 




(25) 
(26) 



Jim V(<p) = 1 



otherwise , 



Pi(<p) = or 



PM = 1 



(27) 



8 ) See contribution to the panel discussion at this meeting by F.H.M. Faisal. 



in which case the constraint is of course not very powerful. In that situation 
it simply means that we have lost too much accuracy in the derivation for that 
particular parameter setting. One should note, however, there is no need to give 
up in that situation since as is evident the expressions for the bounds are by 
no means unique. It should then be quite clear that one can not deduce 9 that the 
bound is useless if one encounters the situation (27). Even more such a conclusion 
seems very much astray in the light of [1,2,4,6], where we presented numerous 
examples for which the bounds are well beyond the values in (27). Sometimes this 
could, however, only be achieved for extremely short pulses. As we pointed out in 
[2] this can be overcome at the cost of having to deal with higher Rydberg states 10 , 
which is a direct consequence of the scaling behaviour outlined in the previous 
section. 

How do typical expressions look like? In [1] we derived for instance the expression 



Vi(ip) = 1 - 
2 



\\(V{x-c(t)e g )-V(x))<p\\dt + 



2E + 6(r) 2 



\\(V(x-c(r)e z )-V(x)M + 



2|&(r)| 
2E + b(r) : 



■\\Pz<p\ 



(28) 



for a lower bound. For given potentials and pulse shapes terms involved in (28) 

may be computed at ease. As stated in [1], it is important to pay attention to 
the fact that (28) is derived for the condition 6(r) 2 /2 > — E = binding energy 11 . 
Such restrictions which at first emerge as technical requirements in the derivations 
usually indicate at some physical implications. In this case it points at the different 
physical situation we encounter when the total momentum transfer is vanishing (see 
also (25)). 

What still needs to be done? Probably it is unrealistic to expect to find a bound 
which is universally applicable and restrictive at the same time, rather one should 
optimize the bounds for particular situations. For instance it would be highly 
desirable to find more powerful bounds for the situations b(r) = 0, c(r) ^ and 
b(r) = c(r) = 0. For the latter case we expect in hindsight from (25) that the 
loss in the estimations may be minimized if in (23) we chose to compare the Stark 
Hamiltonian with the free Hamiltonian —A/2 instead of H — —A/2 + x • E (t) as 
was done in [1]. 



9 ) As was done by J.H. Eberly at the panel discussion at this meeting. 

10 ) This should not lead to the conclusion that bounds in general are exclusively applicable to 
higher Rydberg states, see contribution to the panel discussion at this meeting by M. Gavrila. 
n ) During the panel discussion at this meeting J.H. Eberly exhibited a plot of our result for 
Vl((p) involving a pulse which did not satisfy this condition. As he confirmed to a question from 
the audience his pulse satisfied 6(r) = 0. The conclusions drawn by J.H. Eberly concerning the 
usefulness of bounds based on this plot are therefore meaningless. (See also footnote 9.) 



IMPORTANCE OF PULSE SHAPES 



From our previous discussion it is evident that the physical outcome differs for 
different pulse shapes. However, the fact that a pulse is adiabatically switched 
on or off is not very important, rather the precise values of b{r) and c(r) are the 
determining quantities. In particular the case 

6( r ) = c (r) = (29) 

is very special, since then asymptotically weak stabilization is certain to exist. An 
adiabatically switched on or off pulse sometimes satisfies (29), but this condition 
is by no means identical to it. We found no evidence for stabilization for an adi- 
abatically switched on field when fe(r) ^ 0. To our knowledge the importance of 
(29) was first pointed out by Grobe and Fedorov [11], using intuitive arguments, 
who employed a trapezoidal enveloping function with symmetrical turn-on and 
turn-off time T, which has the nice feature that for T and r being integer cycles 
b(r) = c(r) = and for T half r being integer cycles b(r) = 0, c(r) ^ 0. There- 
after, this observation seems to have been widely ignored in the literature since 
many authors still employ pulses which do not have this property, trading (29) 
for the condition of an adiabatic smooth turn-on or/and turn-off 12 . For instance 
a sine-squared switch on and off with T and r being integer cycles has 6(r) = 0, 
c(r) 7^ 0, an entire sine-squared envelope for r being integer cycles satisfies b{r) = 0, 
c(t) 7^ 0. Using gauBian envelopes or gauBian switch on and no switch off usually 
yields b(r) ^ 0, c(r) 7^ 0. A pulse which has the nice features that it allows a 
theoretical investigation of all possible cases for the values of b{r) and c(r) is the 
triple (5-kick in the form 

f(t) = 5(t) + P ± 5(t - r/2) + f3 2 8{t - t) , (30) 

which we employed in [6]. This pulse obviously satisfies 

6(r) =£ (l+ft + ft/2) and c(r) = E (l + ft/2) (31) 

such that by tuning the constants ft, (3 2 we may realise any desired value of b{r) 
and c(t). 

How do real pulses look like 13 ? The quantity which is experimentally accessible is 
the Fourier transform of the pulse (2) 

POD °° 

E(u) = / E(t)e wt dt = "nw" • (32) 



12 ) As may be supported by numerous publications, this observation appears not to have become 
common knowledge as claimed by M. Gavrila in the introduction to the panel discussion at this 
meeting. 

13 ) We acknowledge that the following argument was initiated, though not agreed upon in this 
form, by an e-mail communication with H.G. Muller. 



with a n being constants. For finite pulses this quantity coincides with the total 
momentum transfer for vanishing frequency uj 

/OO rT 
E(t)dt = / E(t)dt = b(r) . (33) 
-oo JO 

Provided that «o — b(r) = 0, the Fourier transform of the momentum transfer 

/oo 
b(t)e lujt dt (34) 
-oo 

is on the other hand related to the total displacement for vanishing frequency 

/•oo rr 
b(t)dt= / b(t)dt = c(r) (35) 
-oo JO 

such that 

E{uj) = btye™^ - iujb{uj) ~ -iuc(r) + 0{uj 2 ) . (36) 

This means that when the experimental outcome is 

E(uo) = a 2 oo 2 + a 3 uj 3 + a 4 o; 4 + . . . (37) 

the total momentum transfer and the total displacement are zero. Experimentally 
the observed fall off is expected to be even stronger [12]. 



COMPARISON WITH GV-PERTURBATION THEORY 

It is instructive to compare our findings with other standard methods as for 
instance the Gordon- Volkov (GV) perturbation theory. Using now in (23) for H 2 
the Hamiltonian just involving the field and the free particle Hamiltonian in the 
Kramers-Henneberger frame subsequent iteration yields 

U 1 (t,t) = U GV (t,t)-i [ dsU GV (t,s)VU GV (s,t) 

Jt' 

- J ds ds' U GV (t, s') VU GV (s', s)VU GV {s, t) + ... (38) 

where U G y corresponds to the free-particle evolution operator in the KH frame 

U GV (t, t) = e -f(*) e -»(t> e ^)i'. e -<M)^ e -fc(*')p. e ft(f)*e»(*'). (39) 

As was explained in [4] we may use these expressions together with the Riemann- 
Lebesgue theorem in order to obtain the extreme frequency and intensity limit, 
finding (24), (25) and (26). For these arguments to be valid we have to assume 



that the Gordon- Volkov series makes sense, so in particular we have to assume its 
convergence. 

The latter assumption may be made more rigorous when considering the one- 
dimensional delta potential V (x) = —a 5 (x) which is well known to possess only 
one bound state. In that case the problem of computing ionization probabilities is 
reduced to the evaluation of 

V(<p) = l- \Mgv(t)) + (v^(r)>| 2 (40) 

with 

2 poo exp ( — ira 2 ^ — ic(r)ap) 

Mov(t)) = -e-*V dp ) J (41) 

k J-oo (l + (p + 6(r)/a) ) (1+p 2 ) 



= «ShL /j' (s) ™fT' (42) 

Here the only unknown is the function ipj (t) which can be obtained as a solution 
of the Volterra equation 



f°° if e 2 ('- s ) 

(t) = J dp^ GV (p,t) + a^J — J dsip! (s) . (43) 



Iteration of this equation is a well controllable procedure and in [6] we found that 
the series converges for all values of a. The results obtained from the analysis of 
this equation match the results obtained from bounds. 



CONCLUSIONS 

The main outcome of our investigations is that the classical momentum trans- 
fer and displacement caused by a laser pulse on an electron are the essential 
parameters determining the existence of weak asymptotic stabilization. In fact, we 
obtained evidence for stabilization only for pulses for which these two quantities 
vanish at the end of the pulse, i.e., with b(r) = and c(r) = 0. 

Using purely analytical methods, we have shown that, for a wide range of po- 
tentials, namely Kato and one- and three-dimensional delta potentials, we always 
have liniE^oo V(ip) = 1 unless b(r) = and c(r) = 0, in which case the ionization 
probability tends to the lowest order in GV-perturbation theory, which corresponds 
simply to the free particle Green's function (39). Furthermore, for infinite frequen- 
cies, the high-frequency condition of [10] is a way to obtain b(t) = and c(t) = 
for all times. 



Clearly, smooth pulses in general do not necessarily fullfil the above conditions, 
and therefore will not provide a mechanism for stabilization, but just prolong the 
onset of ionization. In fact, we have observed no stabilization for adiabatically 
switched on and off pulses of several shapes, for which analytic expressions for 
lower bounds of ionization probabilities lead to conclusive statements concerning 
the existence or absence of stabilization. 

Therefore, as an overall conclusion: Bounds are useful indeed, also in the 
context of high intensity laser physics! 
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